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The depletion interaction mediated by non-adsorbing polymers promotes condensation and assem¬ 
bly of repulsive colloidal particles into diverse higher-order structures and materials. One example, 
with particularly rich emergent behaviors, is the formation of two-dimensional colloidal membranes 
from a suspension of filamentous fd viruses, which act as rods with effective repulsive interactions, 
and dextran, which acts as a condensing, depletion-inducing agent. Colloidal membranes exhibit 
chiral twist even when the constituent virus mixture lacks macroscopic chirality, change from a cir¬ 
cular shape to a striking starfish shape upon changing the chirality of constituent rods, and partially 
coalesce via domain walls through which the viruses twist by 180°. We formulate an entropically- 
motivated theory that can quantitatively explain these experimental structures and measurements, 
both previously published and newly performed, over a wide range of experimental conditions. Our 
results elucidate how entropy alone, manifested through the viruses as Frank elastic energy and 
through the depletants as an effective surface tension, drives the formation and behavior of these di¬ 
verse structures. Our generalizable principles propose the existence of analogous effects in molecular 
membranes and can be exploited in the design of reconfigurable colloidal structures. 


I. INTRODUCTION 

Suspensions of particles with hard-core repulsive inter¬ 
actions form equilibrium phases that minimize the sys¬ 
tems’ free energy by maximizing their entropy. Since en¬ 
tropy is conventionally associated with disorder, it might 
be expected that hard-particle fluids form structures that 
lack long-range order. However, extensive experimen¬ 
tal work and theoretical models have repeatedly demon¬ 
strated the counterintuitive notion that entropy alone is 
sufficient to stabilize ordered phases of ever-increasing 
complexity. Among other examples, it has been shown 
that entropy can drive formation of 3D bulk crystals in 
suspensions of hard spheres [I], nematic and smectic liq¬ 
uid crystalline phases with hard rods 0 0 , and more 
exotic binary crystals and diverse microphase-separated 
states in mixtures of hard particles mm- 

Recent work has demonstrated that a mixture of 
monodisperse micron-long filamentous bacteriophages 
and non-adsorbing polymers assemble into 2D one- 
rod-length-thick colloidal monolayer membranes mm- 
Colloidal membranes exhibit an exceedingly rich phe¬ 
nomenology. They support a myriad of defects including 
twist domain walls and linear arrays of pores 0 . In- 
creasing chirality induces a transition of flat 2D mem¬ 
branes into ID twisted ribbons, and mixing rods of 
multiple lengths leads to formation of finite-sized col¬ 
loidal rafts that are evocative of similar structures ob¬ 
served in conventional lipid bilayers mm- All of these 
complex mesoscopic behaviors arise from very simple 
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microscopic interactions between constituent particles. 
Filamentous viruses interact only through an effective 
hard-rod repulsion. Similarly, the uncharged dextran 
molecules act as effective Asakura-Oosawa penetrable 
spheres [TlJ □ 3 . From this perspective, the virus par¬ 
ticles and dextran molecules comprise a gas of hard rods 
and hard spheres, and the structures found in colloidal 
membranes must be stabilized by entropic, hard-core in¬ 
teractions m- We formulate a theoretical model based 
purely on such entropic considerations. Our model ex¬ 
plains many known structural features of colloidal mem¬ 
branes and directly relates them to the known entropic 
interactions in rod/polymer mixtures. Furthermore, it 
makes a number of new predictions that are directly ver¬ 
ified by new experimental results. 

Colloidal suspensions are a quintessential model sys¬ 
tem in soft condensed matter physics. They are not only 
interesting in their own right but also provide new in¬ 
sights into the structure and dynamics of diverse phases; 
these insights only depend on the symmetries of the con¬ 
stituent particles and are thus relevant on all length- 
scales. For example, engineering colloidal shapes and 
interactions makes it possible to mimic many processes 
found in atomic and molecular systems, including liquid- 
gas phase separation, wetting, thermal capillary waves, 
crystal nucleation, and the glass transition HUMS. In 
stark contrast to molecular systems, the size of model col¬ 
loids makes it is possible to directly track the positions 
of all the constituent particles, thus yielding important 
information about universal physical processes in vari¬ 
ous condensed matter systems. Conventional fluid mem¬ 
branes, assembled from permanently-linked hydrophobic 
and hydrophilic components, are another interesting and 
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important soft matter system and play an essential role 
in biology m- However, due to our inability to di¬ 
rectly visualize real-time dynamics of lipid bilayers at 
the nanometer scale, many membrane-based processes 
remain poorly understood. Intriguingly, the large-scale 
elastic deformations of colloidal membranes are described 
by the same continuum theories that are used to describe 
conventional lipid bilayers. Based on this observation 
and following the analogy between colloids and molec¬ 
ular substances, we hope that colloidal membranes will 
provide new understanding about universal membrane- 
mediated behaviors. There have been some recent over¬ 
tures in this vein. For example, colloidal membranes per¬ 
mit direct visualization and quantitative characterization 
of liquid raft-like clusters pm, a subject that remains con¬ 
troversial in conventional lipid membranes uni mi- Even¬ 
tual understanding of such complex structures requires a 
theoretical model that relates mesoscopic properties of 
colloidal membranes to the microscopic interactions of 
their constituent building blocks. 

The rest of the paper is organized as follows. In 
Sec. |TT] we briefly review the rich phenomenology of 
colloidal membranes. In Sec. |III[ we introduce a new 
entropy-based theoretical model of colloidal membranes 
and compare our results to known properties of col¬ 
loidal membranes, including static edge fluctuation data 
[Figs. (5jb), [5][d), andj5je)] and twist domain wall re- 
tardance (Fig. [8| [8] OjT Furthermore, we also discuss 
new predictions of our theoretical model, including how 
the structure of the membrane’s edge depends on mem¬ 
brane radius (Fig. [4| and dynamical edge fluctuation 
data [Fig. |5jc)] . These predictions are tested against 
new experimental data. Section |TV| explains the model in 
complete detail and Sec. |V| describes experimental meth¬ 
ods. Finally, we summarize our findings and discuss their 
wider implications in Sec. m 


II. OVERVIEW OF COLLOIDAL MEMBRANES 

Filamentous fd viruses are monodisperse semi-rigid fil¬ 
aments with 880 nm length, 7 nm diameter, and 2.8 pm 
persistence length [21]. When suspended in an aque¬ 
ous solution at increasing concentrations, they undergo a 
transition to an aligned nematic phase characterized by 
long-range orientational order. This isotropic-to-nematic 
phase transition is quantitatively described by Onsager’s 
theory, indicating that viruses repel one another via 
hard-core and electrostatic interactions [21 EH [22]. Fil¬ 
amentous viruses are chiral and form a twisted nematic 
(cholesteric) phase in which the director field rotates with 
a well-defined handedness [23] . For wildtype fd virus, 
the strength of cholesteric interactions is temperature- 
dependent and continuously increases with decreasing 
temperature. A single amino acid substitution in the ma¬ 
jor coat protein leads to the Y21M virus whose cholesteric 
phase has a handedness opposite to that of the wild- 
type m- Mixing wildtype and Y21M viruses produces 


cholesteric phases with intermediate twist pitches; at a 
certain ratio, the mixture exhibits no macroscopic twist. 

The addition of a non-adsorbing polymer, such as dex- 
tran, to a dilute isotropic fd suspension induces virus- 
virus attraction via depletion mi m- The geometry 
of the constituent rods ensures that attractive interac¬ 
tions are strongest for lateral associations, causing the 
viruses to coalesce into one-rod-length-thick, disk-shaped 
mesoscopic clusters [6]. They slowly sediment to the bot¬ 
tom of the glass container, which is coated with a poly¬ 
acrylamide brush penetrable to dextran in order to sup¬ 
press depletion-induced virus-wall attractions m • Over 
a certain range of depletant concentrations, protrusion 
fluctuations induce vertical repulsion between clusters, 
suppressing their face-on association [7.. Consequently, 
such clusters continue to associate laterally, forming large 
equilibrium 2D colloidal membranes that can be millime¬ 
ters in diameter [Fig. [l](b)] . Single molecule tracking 
indicates liquid-like order within a membrane. Twist¬ 
ing of constituent chiral viruses is inherently incompat¬ 
ible with assembly into a layered membrane-like struc¬ 
ture [25]. Consequently, twist can only penetrate into the 
membrane from the edges and is expelled from the bulk. 
Unique properties of the colloidal membrane allow for di¬ 
rect visualization of the twist field and quantitative mea¬ 
surement of the twist penetration length /twist [25 . When 
the membrane radius is much bigger than /twist ^ 1 pm, 
the edge adopts a surface-area-minimizing rounded shape 
with the constituent rods significantly tilting into the 
membrane plane [Fig. 0b)] ; when the membrane radius 
becomes of the order of /twist or smaller, the edge profile 
becomes more square-like and rods do not significantly 
tilt away from the membrane normal [Fig. [ljc)]. Due to 
thermal excitations, membrane edges undergo ripple fluc¬ 
tuations that can be visualized and precisely quantified 
[Fig.0d)j. 

When chirality-inverted Y21M viruses are used instead 
of wildtype fd, rods at the edge twist with the opposite 
handedness, and when the macroscopically achiral mix¬ 
ture of wildtype and Y21M viruses is used, edge-bound 
rods in each membrane have equal probability of twisting 
with one handedness or the other [9]. The achiral mix¬ 
ture exhibits spontaneous symmetry breaking, which has 
been observed in Langmuir-Blodgett films mum , an¬ 
other class of two-dimensional structures with nanoscale 
components, and which has been used in sensors of molec¬ 
ular chirality [28]. Increasing the rod chirality raises the 
free energy of interior untwisted rods while lowering the 
free energy of edge-bound twisted rods, leading to chi¬ 
ral control of edge line tension [9]. At sufficiently high 
chirality, the edge tension approaches zero, and a flat 2D 
disk spontaneously transitions into an array of ID twisted 
ribbons, called a “starfish” [Fig. [lje)]. 

The twist associated with the membranes edge also 
leads to unconventional pathways of membrane coales¬ 
cence [8]. As two membranes of same chirality approach 
each other laterally, the proximal membrane edges can 
partically coalesce and localize 180° of twist to a ID 
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Structures and properties 





FIG. 1. (Color) Overview of colloidal membranes, (a) fd vims particles and dextran molecules act as rod-shaped colloids 
and spherical depletants, respectively, (b) Depleting molecules condense a dilute isotropic virus suspension into a liquid-like 
colloidal monolayer of aligned rods. From left to right, differential interference contrast (DIC) image of circular membranes 
of various sizes, transmission electron microscopy image showing a curved cross-section of the edge of a large membrane, and 
schematic of two large circular membranes of opposite chirality, (c) From left to right, top- and side-view LC-Polscope images of 
a medium-sized membrane, top- and side-view LC-Polscope images of a small membrane, and top- and side-view LC schematics 
of a small membrane. Along with (b), these images illustrate that edges of smaller membranes are more squared, (d) DIC 
images of thermally-excited ripple fluctuations at four different times, (e) DIC images of a temperature induced transition of 
a flat 2D colloidal membrane (left) into a structure with a starfish morphology (right), (f) DIC image (left) and schematic 
(right) of a twist domain wall, or 7r-wall, formed from two partially-coalesced circular membranes, (b) (left), (c), (d), (e), (f) 
Scale bars, 4 pm. (b) (middle) Scale bar, 0.2 pm. 


structure between the membranes; consequently, such 
structures are called it- walls [Fig. [ljf )]. The rods twist by 
180° along the axis connecting the two membranes, from 
one side of the 7r-wall to the other. At the middle of the 
7r-wall, the rods point in the plane of the membranes. 


III. RESULTS 

A. Circular membranes 

In our model, we treat the membrane as a continuous 
fluid composed of rods at constant density. Once the 
membrane is stably formed, we assume it does not ex¬ 
change rods with the surrounding solution; thus, its vol- 
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lime is fixed. The membrane structure is characterized 
by two coarse-grained degrees of freedom available to the 
rods: a twist angle 0(x) about an axis in the membrane 
plane and a root-mean-square amplitude b(x) of height 
fluctuations perpendicular to the membrane plane. Per¬ 
pendicular fluctuations increase the effective thickness of 
the membrane, and instead of using 6(x) directly, we will 
develop a microscopic theoretical model and present its 
results using the coarse-grained membrane half-thickness 
h(x) = £cos0(x) + b(x), where t is the half-length of the 
virus. 

The model free energy is comprised of three entropic 
components. The first term is the Frank free energy 
that disfavors bend elastic distortions of rods within a 
membrane while favoring local twisting of rods at their 
naturally-preferred wavenumber q m-, it depends pre¬ 
dominantly on 9. All experimental results are obtained 
using wildtype virus suspensions, which favor left-handed 
twist [2T] . The second term is associated with free vol¬ 
ume accessible to the depleting polymer due to the pres¬ 
ence of the membrane; it depends on the thickness pro¬ 
file h [see Fig. |2ja)]. Aside from the constant volume of 
the incompressible membrane, the excluded volume is ap¬ 
proximately its surface area times the depletant radius; 
thus, this term acts as an effective surface tension energy. 
Its magnitude is proportional to the depletant concentra¬ 
tion and to the temperature. The third and final term 
accounts for the entropy associated with rods protrud¬ 
ing from membranes into the surrounding volume occu¬ 
pied by the depleting polymer, a phenomenon reported 
in Ref. [7]. Protrusion of each rod increases the effective 
surface area of the membrane, which decreases the vol¬ 
ume accessible to the depletant molecules. The preferred 
magnitude of rod height fluctuations bo is determined by 
a trade-off between rod entropy, which prefers larger 6, 
and the depletion effect, which tends to minimize b [see 
Fig. [2^b)] . In our system, the preferred magnitude of this 
effective surface roughness is very small —bo t —but 

the energetic cost of deviations from this value depends 
on the rod angle 9 [see Fig. [2jc)] . When 9 « 0, rods 
are packed more closely in the plane of the membrane, 
assuming a constant perpendicular distance £ between 
rods. Thus, rod fluctuations produce surface roughness 
on a smaller length scale, which creates more effective 
surface area and costs more energy. In this case, b = bo is 
strongly preferred, so h ~ t cos 9 and rod entropy can be 
ignored. When 9 ~ 1, rods are spaced farther apart in the 
plane of the membrane, leading to fluctuation-produced 
surface roughness on larger length scales. These longer- 
wavelength fluctuations resist b = bo more weakly, so h 
may differ significantly from t cos 9. In a similar fashion, 
manipulating the surface roughness of larger colloids can 
tune their depletion-induced interaction [30] [31]. In sum¬ 
mary, the rod fluctuation term couples h to t cos 9 with a 
0-dependent coupling strength. To obtain the membrane 
structure, we minimize the total free energy over 0(x) 
and /i(x). At the center of the membrane, the membrane 
is fixed to be one-virus-length thick, while there are no 
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FIG. 2. (Color) Depletant contributions to the membrane free 
energy. Depleting molecules (dextran polymers in our system) 
are illustrated in green, and the volume excluded to them by 
the membrane is illustrated in blue, (a) The volume that is 
excluded to the depleting polymer due to a a smooth mem¬ 
brane is comprised of the volume of the membrane itself and 
to a first order the membrane surface area times the depletant 
radius. h(x) is the coarse-grained half-thickness of the mem¬ 
brane. (b) Rod height fluctuations produce an effective sur¬ 
face roughness that increases the effective surface area of the 
membrane and contributes extra excluded volume: (top) an 
idealized configuration with no height fluctuations and (bot¬ 
tom) an idealized configuration with small-wavelength height 
fluctuations of amplitude b — h — t cos 9 that increase the ex¬ 
cluded volume. The dotted line represents the excluded vol¬ 
ume of the configuration without height fluctuations, (c) The 
free energy density of rod height fluctuations are calculated 
in the mean-field limit by considering a single protruding rod 
amid a membrane of constant local thickness: (top) rods at 
small tilt angle 9 and (bottom) rods at large 9. The mag¬ 
nitude of this free energy density decreases with increasing 
9 because tilted rods are less dense in the membrane plane 
by a factor of cos0, assuming a constant perpendicular dis¬ 
tance £ between rods. In other words, the surface roughness 
lengthscale in the ^-direction is proportional to l/cos0, or 
equivalently, the extra effective surface area created by fluc¬ 
tuations is proportional to cos 9. Theoretical details are given 
in Sec. UVAl 


height constraints at the membrane edge. 

In order to obtain quantitatively meaningful results, 
we use parameter values that are extracted from rele¬ 
vant experimental measurements when possible (Table]!]). 
Five parameters, whose values are neither experimen¬ 
tally controlled nor directly measured, are allowed to 
vary as fit parameters: the characteristic depletant size 
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TABLE I. Membrane parameters and their values. 


Parameter 

Variable 

Experimental value 

Reference (s) 

Theoretical fit value 

Vims half-length 

t 

440 nm 

ED 

same 

Temperature 

T 

0-60 °C 

experimental 

same 

Depletant concentration 

n 

35-51 mg mL -1 

experimental 

same 

Depletant radius 

a 

~25 nm 


31 nm 

Nearest-neighbor vims distance 

£ 

12 nm 

unpublished 15 

same 

Frank elastic constant 

K 

0.5 pN 

pyqJcTI '—' 

2.8 pN 

Preferred twist wavenumber 

q(T) 

0.5 pm -1 — T/60 °C 

[90 

2.5 pm - Vl-T/120°C 

Vims birefringence 

An 

0.0087 ±0.0007 

25 o 

0.0065 


a Hydrodynamic radii for dilute solutions of 500 kDa dextran, whereas our experiments are in the semidilute regime. 
b Unpublished data extracted from X-ray scattering. 

c Measured in the bulk cholesteric phase with fd virus concentration lOOmgmL -1 , which is lower than the membrane virus 
concentration 230mgmL —1 estimated from the experimentally-measured nearest-neighbor virus distance £. 
d Assuming that the nematic order parameter in membrane is 1. Membrane virus concentration 230mgmL —1 estimated from the 
experimentally-measured nearest-neighbor virus distance £. 


a, the Frank elastic constant K in the one-constant ap¬ 
proximation, the amplitude and transition temperature 
of the temperature-dependent twist wavenumber g(T), 
and the virus birefringence An. In our theory, we main¬ 
tain the experimentally-measured square-root behavior 
of q{T) (see Supplementary Material of [9]). The Frank 
elastic constant can be written in dimensionless form as 
k(T) = K/natT , a ratio between the influence of Frank 
elasticity and that of depletion. Presumably, K depends 
on temperature in a complicated fashion, as measured 
for a variety of lyotropic and thermotropic liquid crys¬ 
tals [35H38] . but we ignore this effect. 

We first use our theoretical model to determine how 
membrane structure depends on its radius. We use 
cylindrical coordinates and assume circular symmetry 
[Figs. §a) and[3fb)]. For convenience, we use the re¬ 
verse radial coordinate Ar, which originates at the mem¬ 
brane’s edge and takes positive values towards the center 
of the membrane. 0 is the twist angle about the local ra¬ 
dial axis. Figure [3f c) plots the vertical membrane profile 
for membranes with very large radii and varying Frank- 
to-depletion ratios k and twist wavenumbers q. For all 
conditions, h ~ tcosO , indicating that 0 is sufficiently 
small to suppress rod height fluctuations. Thus, rod en¬ 
tropy does not contribute significantly to the structure of 
the membrane’s edge. First, consider the q = 0 profiles 
in Fig. |3fc) corresponding to a macroscopically achiral 
rod mixture. When k is greater than a critical value 
k c = 1 , the untwisted configuration with 6 = 0 is fa¬ 
vored. When k < fc c , depletion drives spontaneous chiral 
symmetry breaking into a twisted configuration with ei¬ 
ther 6 > 0 or 0 < 0. In the k —> 0 limit where only 
depletion exists, the vertical edge profile is semicircular 
to minimize the membrane surface area. Now, consider 
the q = 2.2 pm -1 case in Fig. |3fc) corresponding to a 
chiral rod mixture. Twisted configurations of one hand¬ 
edness (here, 6 > 0 for q > 0) become favored at all 
k. In the depletion-dominated regime k <C 1, the ver¬ 
tical edge profile again approaches a semicircle. In the 
Frank-elasticity-dominated regime k 1, the rod twist 


decays with penetration length /twist ~ Vk/t , in analogy 
to the way that twist penetrates into a smectic phase. 
Calculations of k c and / tw ist are provided in Appendix |A| 

In addition to describing edges of large membranes, 
our theoretical model also describes how edge profile 
varies with decreasing membrane diameter. To test these 
predictions, we measure the retardance of different-sized 
membranes using quantitative LC-PolScope microscopy, 
which directly reveals the twisting of rods away from the 
membrane normal. When polarized light passes through 
a birefringent material, the components corresponding to 
the dielectric tensor eigenvectors—the ordinary and ex- 
traodinary waves—propagate at different speeds. The re¬ 
sulting phase difference between these components mul¬ 
tiplied by the wavelength of the light is the retardance 
D. For a uniaxial crystal of constant thickness, retar¬ 
dance can be calculated as D = A nh sin 2 9 [39], where 
An is the birefingence. For membranes of various radii, 
we calculate D(Ar) with the fit values given in Table [I] 
and the approximation h = tcosO , since our results in 
Fig-! c) demonstrate that rod fluctuations b are insignif¬ 
icant for membrane edges. We use the same parameter 
values for all membrane sizes; only the radius changes. 
The radially-averaged edge retardance profiles measured 
for membranes of various radii match well with our theo¬ 
retical predictions [Fig.[4fb)]. These results demonstrate 
that rods are less tilted at the edges of smaller mem¬ 
branes compared to those of larger membranes [insets of 
Fig |4](b)] , consistent with observations that larger mem¬ 
branes appear on side-view to have rounded edges while 
smaller membranes have squared-off edges [Figs. Ilf b) and 

0c)]. 

With detailed understanding of the membrane’s edge 
structure, we next study its fluctuations, which are 
clearly visible and easily quantified with optical mi¬ 
croscopy [Fig. [ljd)]. In the large membrane limit, we 
ignore curvature of the edge and, with Cartesian coor¬ 
dinates, place the very edge at x = 0 [Fig. 0 a)]. 0 is 
now the twist angle about the x-axis. Using the previ- 
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FIG. 3. (Color) Vertical edge profile of a large membrane and its dependence on Frank elasticity and chirality, (a) Perspective 
and (b) cross-section schematics show parametrization of the membrane edge profile and the cylindrical coordinate system, 
(b) shows rods that intersect the light blue plane in (a). A r is a reverse radial coordinate where A r — 0 corresponds to the 
membrane edge, h is the membrane half-thickness and 0 is the rod tilt angle, t is the half-length of the rods, (c) Calculated 
edge profiles of a large membrane (radius R^t) with various Frank-to-depletion ratios k from left to right and preferred twist 
wavenumbers q from top to bottom. In all cases, h (blue) is almost indistinguishable from t cos# (cos 0 in red). Note that for 
q — 0 and k < 1, 0 ^ 0, demonstrating spontaneous symmetry breaking into a configuration with one handedness (6 > 0) or 
the other (0 < 0). For q — 0 and k > 1, the untwisted state with 0 = 0 has lowest energy. Experimental conditions listed in 
Table [T] are closest to k — 0.85 and q — 2.2 pm -1 . 







A r [pm] 


FIG. 4. (Color) Membrane edge retardance. (a) 2D LC-PolScope birefringence map of a large circular membrane with retardance 
represented as pixel brightness. The dotted green line approximately corresponds to the range of Ar’s plotted in (b). Scale bar, 
4 pm. (b) Retardance values D for circular membranes of various sizes. The points indicate experimental data at temperature 
T = 22 °C and depletant concentration n — 45mgmL _1 . The lines indicate theoretical results calculated with these parameter 
values and those described in Table [T] giving k = 0.85 and q = 2.2 pm -1 . Membrane radii range from 5.1pm (top left) to 
0.45pm (lower right). The insets show the calculated membrane profile /i(r), plotted with an aspect ratio of 1. Tick marks 
signify 1pm increments. tcosO , not shown, is strongly coupled to h in all cases. 


ously discussed model, we first calculate h{x) and 0{pc) 
for a flat edge. We then introduce a small ripple at the 
edge with the tangent angle a(y) that perturbs the rod 
configuration as in Fig. [5|a). The unperturbed configu¬ 
ration along lines parallel to the y -axis is mapped onto 
curves with the same tangent angle a(y), and the rod ro¬ 
tation axis for 0 is always perpendicular to these curves. 


See Section |IVB| for a mathematical description of this 
ripple ansatz. We can write a(y) in terms of Fourier com¬ 
ponents a pi where p is the ripple wavenumber. To lowest 
order in these Fourier components, the relative free en- 
ergy per unit length is f = \ E P (7[M] + 0}p 2 )a 2 p , 

where the line tension 7 [ft, 0\ and the edge bending mod¬ 
ulus K,[h, 0] are functionals of the flat edge configuration. 
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FIG. 5. (Color) Line tension analysis of the membrane edge. 
The points indicate experimental data at various temper¬ 
atures T and depletant concentrations n. The lines indi¬ 
cate theoretical results calculated for corresponding param¬ 
eter values and those described in Table [T] (a) Schematic of 
the membrane ripple ansatz through which line tension and 
line bending modulus are calculated. a(y) is the angle be¬ 
tween the ripple tangent vector and the y-axis. (b) Thermal 
fluctuation amplitudes (a: 2 ) and (c) autocorrelation decay 
timescales 1 /uj p of ripple fluctuations for depleting concen¬ 
tration n 45mgmL _1 and various temperatures T. The 
theoretical plots of 1 /cj p use the fit value for the ID mem¬ 
brane edge viscosity rjin = 300mPaspm 2 . (d) Line tension 
7 and (e) its relative temperature-dependent behavior as a 
function of temperature for various n. For each n, y 0 ° c is the 
line tension extrapolated to T = 0 °C. 


The line tension describes the energetic cost of having 
an edge interface, and the edge bending modulus arises 
from the rod director’s bend distortion introduced by the 
ripple. By equipartition and viscous hydrodynamics, we 
obtain the fluctuation spectra 


fc bT 1 _ rj ip 

7 + Hip 2 ’ UJp 7 + Hip 2 


(<a 2 ) is the average fluctuation amplitude of Fourier mode 
p and uj p is the temporal autocorrelation decay constant 
of Fourier mode p, as found in the temporal autocorre¬ 
lation function (a p (t)a p ( 0)) = {a^)e~ Upt . The decay of 
fluctuation correlations arises from dissipative forces, the 
most significant of which are membrane viscous stresses 


since we expect the membrane to be much more viscous 
than the solvent, pm is the one-dimensional (ID) viscos¬ 
ity of the membrane edge. 

Using the fit values in Table [T] describing the membrane 
edge, our theoretical model predicts values for 7 and k, 
which determine the fluctuation spectra (<a 2 ) and 1 /u) p . 
These predictions can be tested experimentally, and the 
value of 7 can be extracted from the low-p limit of (<a 2 ). 
The experimental and theoretical spectra match well over 
a variety of temperatures [Figs. [5jb) and[5][c)]. These cal¬ 
culations still assume h = £cos$, since Fig. p3[c) demon¬ 
strates that rod fluctuations b are insignificant for mem¬ 
brane edges. The ratio between 1/lj p and (<a 2 ) appears 
constant for all measured values of p —in agreement with 
our theory—and gives a value for 7710 ~ 300 mPa s pm 2 . 

We expect the 3D membrane viscosity 77 to be strongly 
inhomogeneous and anisotropic at the edge due to the 
large aspect ratio of the rods. For instance, during a rip¬ 
ple fluctuation, rods oriented more vertically may slide 
past each other more easily than those tilted more hor¬ 
izontally. To roughly estimate the magnitude of 77, we 
write 771 D 00 J dxdzp ~ Aip where A ~ 2£/ twist ^ 2t 2 is 
an estimated cross-sectional area of the membrane edge 
participating in these ripple fluctuations. As calculated 
in Appendix [aJ /twist ~ Vkt is the twist penetration 
depth, and the parameter values provided in Table [I] sat¬ 
isfy k ~ 1. This gives 77 ^ 800mPas, much greater than 
the solvent viscosity, which is rj s ~ 3mPas for 5 w% 
500 kDa aqueous dextran uni. 

Measurements and calculations of the line tension 7 
show good quantitative agreement over a variety of tem¬ 
peratures T and depletant concentrations n [Fig. |5jd)]. 
For all 72 , 7 decreases as T is reduced. If we measure 
7 relative to its value at a standard temperature, say 
T = 0°C, the line tensions for different n all collapse 
onto a single curve [Fig.j5je)], indicating that the relative 
effect of temperature change on 7 is independent of the 
depletant concentration. In Sec. |IVB[ we see how these 
effects arise naturally in our model via a ^-dependent chi¬ 
ral term in the line tension. Colloidal membranes assem¬ 
bled from chiral rods are inherently frustrated, because 
the particles cannot simultaneously twist locally and as¬ 
semble into a monolayer globally. Consequently, twist 
is expelled from the membrane interior and localized to 
its edges. Note that q(T) is a monotonically decreasing 
function of T (Table [T|). Decreasing the temperature in¬ 
creases q(T) and lowers the free energy of edge-bound 
twisted rods, leading to chiral control of edge line ten¬ 
sion [9]. 


B. Starfish morphological transition and 7r-wall 
structure 

We now apply our theory to explain more exotic 
structures found in colloidal membranes. For example, 
when circular membranes are subjected to a tempera¬ 
ture quench, the line tension decreases significantly and 
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FIG. 6. (Color) Starfish instability transition for large mem¬ 
branes. (a) Starfish arms grow from unstable ripple fluctu¬ 
ations in our theory. For large membranes, we can take the 
Cartesian limit and ignore the curvature of the edge, (b) 
The shaded region indicates unstable ripple wavenumbers p 
calculated for preferred twist wavenumbers q and constant 
Frank-to-depletion ratio k = 0.85. We take p to be contin¬ 
uous, corresponding to the infinite membrane size limit; for 
finite-sized circular membranes, continuity permits only cer¬ 
tain values of p, namely multiples of the inverse radius. The 
inset plots the same results on a log-log scale to demonstrate 
that as q increases past a critical q c ~ 3 pm -1 , the range of 
unstable p’s grows as a power law with exponent 1/2. 


fluctuations at the edge increase in amplitude. For suf¬ 
ficiently low T, the circular membrane becomes unstable 
and grows arms of twisted ribbons along its entire periph¬ 
ery [Fig. [lje)]. In our model, these starfish arms arise 
from the aforementioned ripple fluctuations [Fig. [ 6 ][a)]. 
As the temperature decreases, the chiral wavenumber 
q(T) increases and lowers the line tension 7 . For suffi¬ 
ciently large g, 7 becomes negative and long-wavelength 
ripple modes along the membrane circumference become 
unstable, which presumably grow and twist into starfish 
arms. Figure ^b) plots the range of unstable wavenum¬ 
bers p, measured around the circumference, as a function 
of chiral wavenumber q for constant k = 0.85. Above a 
critical q c ~ 3 pm -1 , low-p modes become unstable. An 
instability with p « 1 pm -1 in a membrane of radius 


R « 5 pm corresponds to a five-armed starfish struc¬ 
ture as depicted in Figs. [lje) and | 6 ][a), so the order of 
magnitude of unstable p’s calculated in Fig. ^b) follows 
expectations. Note that changing the temperature also 
changes fc(T), but the effect is qualitatively insignificant. 
The transition from a circular membrane to a starfish 
structure is reversible, so reheating to a positive 7 drives 
the edge-length-maximizing starfish structure to decrease 
its edge length and become circular again [9 . 


We also use our theoretical model to quantita¬ 
tively explain another prominent and experimentally- 
characterized feature of colloidal membranes: the 7r-wall. 
Observations on the assembly pathways and structure of 
7 r-walls were described in Sec. [IT] and Fig. Ba Briefly, 
two membranes of the same handedness can partially co¬ 
alesce into a single membrane and trap a twist domain 
wall, or 7r-wall, through which the rod director twists by 
180°. To investigate these structures theoretically, we 
use Cartesian coordinates [Figs. B a ) and B b )] and set 
6(x = 0) = 7 r /2 so the rods at the middle of the 7 r-wall 
lie completely in the membrane plane. Fig.[7|[c) plots the 
thickness profile over a range of Frank-to-depletion ratios 
k and chiral twist wavenumbers q. In all cases, h is much 
greater than tcosO at the middle of the wall, since the 
coupling that sets h « tcosO becomes very weak when 
0 « 7 r/ 2 . Remember, h = t cos $ + 6 , where b is the ampli¬ 
tude of rod fluctuations perpendicular to the membrane. 
This means these rods undergo position fluctuations in 
the ^-direction that are many times larger than both their 
projected height tcosO and their diameter d « 0 . 02 1. 
Such a phenomenon would require rods to pass through 
each other, which is theoretically allowed because we ig¬ 
nore rod-rod interactions, but we wish to interpret this 
result physically, h t cos 0 and h^> d indicate that the 
membrane is thicker than multiple layers of tilted rods, 
so these large fluctuations may be physically manifested 
as rods stacking on top of each other. The addition of 
repulsive rod-rod interactions may further increase the 
thickness of the 7 r-wall. As for the 7 r-wall profile, sim¬ 
ilarly to the membrane edge, the depletion-dominated 
regime fc«l leads to a circular profile and the Frank- 
elasticity-dominated regime k 1 leads to slow rod 
twist decay [Fig. [7|c)]. q does not significantly affect 
the 7r-wall profile among the parameter values explored; 
indeed, it appears in a ^-dependent free energy term that 
can almost be integrated to the bound ary, an d 0(x = 0 ) 
and Q(x -T 00) are fixed (see Section IVA). With the 
same parameters used to calculate the edge structure 
(Table 0, the calculated optical retardance of the 7r-wall 
quantitatively matches the experimentally measured pro¬ 
file (Fig. h . 
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FIG. 7. (Color) Vertical 7r-wall profiles and their dependence on Frank elasticity and chirality, (a) Perspective and (b) cross- 
section schematics showing parametrization of 7r-wall profile and Cartesian coordinate system, (b) shows rods that intersect 
the light blue plane in (a), h is the membrane half-thickness and 0 is the rod tilt angle, t is the half-length of the rods, (c) 
Calculated vertical 7r-wall profiles for various Frank-to-depletion ratios k from left to right and preferred twist wavenumbers q 
from top to bottom. In all cases, h (blue) is almost indistinguishable from tcosO (cos 0 in red) away from x = 0. Near x — 0, 
h approaches a finite mid-wall value while cos# approaches 0. Insets highlight the profile near x = 0. Experimental conditions 
listed in Table [i] are closest to k = 0.85 and q = 2.2 pm -1 . 


IV. THEORETICAL DEVELOPMENT 


A. Membrane parametrization and free energy 


We treat the membrane as a continuous medium com¬ 
posed of rods at constant density, and we fix the num¬ 
ber of rods in the membrane by fixing the membrane 
volume. The coarse-grained rod twist angle 0(x), rod 
height fluctuation amplitude b(x), and membrane half¬ 
thickness are related by ft(x) = t cos#(x) + 6(x), where t 
is the half-length of the virus. We will first develop the 
model assuming a circularly-symmetric membrane of ra¬ 
dius R and using cylindrical coordinates in which /i(r), 
b(r), and 0(r) only depend on the radial coordinate. 

We model the rods as liquid crystals whose orientations 
are described by a Frank elastic free energy [29]. In a cir¬ 
cular geometry, the rods point in the ^-direction but can 
twist with angle 0 in the azimuthal direction [Figs. la) 
and^b)]. Using the one-constant approximation, the 
free energy is: 


ftv.nk = K j d 2 xfi[(V ■ n) 2 + (V x n)= 


/ 


= 27 tK / dr h 


— 2gn • V x n] 
r(d r 0) 2 + sin 20 d r 0 - 


— 2 qrd r 0 — q sin 20 


(i) 


sin 2 0 


( 2 ) 


K is the 3D Frank elastic constant and q is the preferred 
twist wavenumber associated with intrinsic chirality of 
the constituent rods, n(r) = sin 0{r)(j) + cos#(r)z is the 


nematic director. The q term breaks chiral symmetry, 
such that for q > 0, twisted membranes with d r 0 > 0 
have lower energy than those with d r 0 < 0. When q = 
0, the total free energy is invariant under the chirality 
inversion 0^—0. 

The depletant polymers act to minimize the volume 
excluded to them by the membrane. For polymers small 
compared to the dimensions of the membrane, this ex¬ 
cluded volume is approximately Vo + a A, where Vo is 
the volume of the membrane, A is the surface area of 
the membrane, and a is the characteristic depletant ra¬ 
dius |4l] [see Fig. [2j[a)] . Vo is constant, so depletion serves 
as an effective surface tension. Consequently, the free en¬ 
ergy is given by: 


Udep — 2Tiak^T 


- Aimak bT 


J d 2 xi/l + (V/i) 2 + J dl h 

" r R , _ 

/ drry/l + (d r h) 2 + Rh(R) 

Jo 


(3) 


,(4) 


where n is the depletant concentration, T is the temper¬ 
ature, and ks is the Boltzmann constant. J dl indicates 
an integral over the membrane edge boundary. 

Finally, we allow rods to fluctuate perpendicularly to 
the membrane plane. In general, these fluctuations have 
complicated, non-linear effects on the free energy, but for 
simplicity, we only consider fluctuations of single rods 
and ignore their interactions and correlations [7]. When 
a single rod at small tilt angle 0 protrudes by a small per¬ 
pendicular distance 2 above a flat coarse-grained mem¬ 
brane surface, it introduces an additional spherical cap of 
volume 7T az 2 that is excluded to the depleting polymers 
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FIG. 8. (Color) 7r-wall retardance. (a) 2D LC-PolScope bire¬ 
fringence map of two circular membranes joined through a 
7r-wall with retardance represented as pixel brightness. The 
dotted green line approximately corresponds to the range of 
rr’s plotted in (b). Scale bar, 4 pm. (b) Retardance values 
D. The points indicate experimental data at temperature 
T = 22 °C and depletant concentration n — 45mgmL _1 . The 
lines indicate theoretical results calculated with these param¬ 
eter values and those described in Tabled] 


[see Fig. [2jc)] . Meanwhile, these protrusions are entrop- 
ically favored by the rods. For a distribution of vertical 
rod displacements p(z), the fluctuation free energy for a 
single rod is a sum of depletant and rod entropy contri¬ 
butions: 

Single = 7rnak B T j d zp(z)z 2 + k B T J dzp(z) log p(z). 

Minimizing this free energy yields Po(z) = 
(27t6o) - 1 / 2 exp(— z 2 /26q), where bo = (27raa) -1 / 2 . 

If all rods were to fluctuate with the preferred am¬ 
plitude bo, then the membrane half-thickness h and rod 
angle 6 would be exactly related as h = t cos 0 + bo . How¬ 
ever, in certain structures such as the mid-planes of i r- 
walls, the Frank and depletion free energies favor pro¬ 
files h(x) and 0{pc) that significantly deviate from this 
relationship. To propertly describe these structures and 
account for the energetic cost of h ^ tcosO + bo, we 
calculate the free energy of Gaussian rod fluctuations 
of amplitude b ^ bo- Using the distribution p(z) = 
(27t 6 2 ) -1 / 2 exp(— z 2 /2b 2 ), the single-rod free energy be¬ 
comes F^ngie = 27 Tnak B T(b — bo) 2 to leading order in 
b — bo- To coarse-grain this expression, we multiply by 
the rod density and integrate over the membrane area. 
For simplicity, we assume the rods are packed hexago- 
nally and maintain a constant perpendicular distance £ 
between nearest-neighbors. In the small 0 limit, the area 


occupied by each rod is \/3£ 2 / cos 0. Our final expression 
for the rod fluctuation free energy is 

8ir 2 nak B T f R 2 

r rod = - 7 = - / drr cos# [h — tcosO — bo) , (5) 

v3 £ 2 Jo 

where we have written b in terms of h and 0. This term 
allows h to deviate from t cos 6+bo with an energy penalty 
corresponding to the magnitude of the deviation. Heuris- 
tically, the energy penalty is proportional to cos 0 because 
at higher 0, the rods are spaced farther apart in the plane 
of the membrane, so height fluctuations of individual rods 
induce less roughness at the membrane surface [Fig. [2j[c)] . 

We minimize the total free energy with volume- 
conserving Lagrange multiplier A 


R — T Fp ran k T Rrod T A 



dr rh 


(6) 


over h(r) and 0(r) to obtain the edge profile. The bound¬ 
ary conditions are ft(O) = t + bo and 0(0) = 0; h(R) and 
0(R) are free. 

Equation [ 6 ] simplifies for large membranes when R is 
much greater than the penetration depth of edge twist 
/twist; the edge becomes essentially straight. We can then 
study the profile of a twisted membrane formed from 
an untwisted rectangular membrane of length L y —* oc 
along the ^/-direction and length 2 L x <C L y along the x- 
direction. We allow the membrane profile to vary along 
the ^-direction and impose reflection symmetry about 
the midline x = L x where the rods are perpendicular 
to the membrane (analogous to r = 0 for the original 
circular geometry). We are interested in the edge pro¬ 
file at x = 0. In this setup, each free energy integral 
becomes its Cartesian version, with Fp r ank losing bend 
distortion terms that arise from a circular geometry. In¬ 
stead of a Lagrange multiplier term, however, volume 
conservation can be directly enforced in the following 
way. The volume of the half of the untwisted membrane 
between x = 0 and x = L x is Vo = 2 (t + bo)L x L y . The 
change in volume brought about by a varying h(x) is 
AU = 2 L y Jq X dx [h(x) — (t + bo)]. To compensate for 
the lost volume, we introduce extra volume at the mem¬ 
brane midline where h(x) = t + bo by adding a width 
A L x of untwisted rods; volume conservation requires 

A L x = -AV/[2(t + b 0 )Ly\ = f 0 Lx dx [1 - h{x)/(t + b 0 )j. 

This extra width increases the half-membrane’s surface 
area by AH = 2L y AL x and, since depletion free energy 
is proportional to surface area, contributes the additional 
term naksTAA to F^ ep . Ignoring a constant term pro¬ 
portional to L x L y , the total free energy becomes 

2 ^ tT ,=1 ' dl ( V,1 + (W "^) + '‘ (0 > 

kt f Lx 

+ — dxh [(< d x 9 ) 2 + 2 qd x 9] 

z Jo 

27r r-^x 

H——— / dx cos 0{h — t cos 0 — bo)\j) 
v3 £ 2 Jo 
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where again, k = K/naiksT. Strictly speaking, the in¬ 
tegrals in the last two terms should extend from 0 to 
L x + A L x , but the contributions to the integrals from L x 
to L x + A L x are zero because d x 0 = 0, 6 = 0, and h = 
t + bo in the interior of the membrane. Comparing Eqs.|6] 
and [7| the additional surface area term is analogous to 
a Lagrange multiplier with value naksT/(t + bo), the ef¬ 
fective osmotic pressure exerted by the depletants on the 
membrane. Also, since this Cartesian parametrization 
implicitly inverts the membrane orientation compared to 
the cylindrical parametrization (instead of decreasing r, 
increasing x moves into the interior of membrane), the 
g-term in Eq. [7] has the opposite sign of the g-terms in 

Eq-H 

For membrane edges calculated in Fig. |3jc), httt cos 0 , 
which means rod height fluctuations b are strongly sup¬ 
pressed. This motivates simplification of the free energy 
by taking the infinite coupling limit in which F ro d en¬ 
forces h = tcosO + bo and therefore disappears from 
the free energy. Using values in Table [IJ we calculate 
bo ~ 0.03£ and make the further approximation that 
these protrusion fluctuations contribute only a small frac¬ 
tion to the membrane thickness and can thus be ne¬ 
glected: bo = 0. Numerical calculations of all edge prop¬ 
erties fixing h — t cos 6 are indistinguishable from those 
using the full theory. Thus, the precise form of F ro d, 
whose derivation required many assumptions, does not 
matter for membrane edges as long as it strongly cou¬ 
ples h to tcosO. This simplification permits derivation 
of some analytical results, including an investigation into 
spontaneous chiral symmetry breaking for q = 0, which 
are given in Appendix |A| 

For 7r-walls, we use Eq. [7] without the boundary de¬ 
pletion term proportional to h(0 ) because x = 0 is the 
middle of the wall and no longer an edge boundary. The 
rods there must lie in the membrane plane, so we gain 
the extra boundary condition 0(0) = 7r/2. Now 0 is fixed 
at both boundaries, so if h were enforced to be a func¬ 
tion of 0 like h = tcosO, the g-term could be integrated 
to a constant and the profiles would not depend on g. 
However, unlike their counterparts at edges, h and 0 are 
independent near x = 0, where calculations show that 
the vertical mid-wall profile satisfies h tcosO', thus, 
the membrane structures depend slightly on q [Fig.^c)]. 
This independence arises due to the angle-dependent cou¬ 
pling strength of F ro d, which has a factor of cos0 in the 
integrand (Eq. [5|. Away from the middle of the wall, 
cos 0^1 and deviations from h = t cos 0 + bo are costly 
for F roc j. As cos0 approaches 0, these deviations cost 
less energy in F ro ^, so other terms such as F ^ ep (Eq. [4] 
without the boundary term) gain influence on the pro¬ 
file configuration. The competition between F ro d, which 
prefers h to decrease with cos0 towards the middle of 
the wall, and F^ ep , which prefers a constant h , sets the 
mid-wall thickness. 

It is worthwhile at this point to compare our theory 
with an alternative one, which we will refer to as the 
KM theory after its developers Kaplan and Meyer Elsa, 


that also produces results in very good agreement with 
experimental observations. First, it should be empha¬ 
sized that the philosophical approaches of the two the¬ 
ories are different. Ours can be viewed as a minimalist 
theory based directly on entropic interactions induced by 
dextran depletants and to a lesser extent by the viruses 
themselves. The KM theory, in the grand tradition of 
liquid-crystal physics, is phenomenological at its core. It 
introduces an order parameter T, inspired by that de¬ 
scribing order in a 3D smectic, that describes the transi¬ 
tion from rods oriented predominantly perpendicular to 
the membrane plane (“smectic” phase with f / 0 ) to 
rods oriented predominantly parallel to the plane of the 
membrane (“cholesteric” phase with T = 0). Though 
the introduction of T provides a useful and predictive 
theory, it is not clear how it could be measured. The 
KM theory also introduces terms in the free energy that 
are not directly present in our theory: one measuring 
the energy cost of surface curvature and two providing 
a favored relative orientation of the surface normal N 
and director n at the top and bottom membrane sur¬ 
faces. However, the term proportional to —h/(t + bo) 
in the Eq. [7] version of our theory provides a preference 
of 0 = 0 , i.e., the director prefers to be parallel to the 
layer normal. More generally, the Lagrange multiplier 
term in Eq. [ 6 ] provides this preference. Naturally, the 
KM theory has more free parameters than the five of 
our theory: depletant size, Frank elastic constant, twist 
wavenumber amplitude and transition temperature, and 
virus birefringence (Table [I]). In spite of these differences 
between the two theories, they share some common fea¬ 
tures: They both employ the Frank free energy with a 
term favoring twist to describe the energetics of director 
deformations, and they both introduce a term favoring 
h = tcosO (when bo can be ignored in our theory) with 
a coefficient (cos 0 in our case and |T | 2 in the KM case) 
that vanishes at a 7r-wall when 0 = 7 r/ 2 , importantly al¬ 
lowing h to differ from t cos 0 with no direct energy cost 
at that point. 

KM pursues a different approach to boundary condi¬ 
tions than we do. They impose the condition 0(R) = n/2 
at the free edge of a circular membrane, whereas we al¬ 
low the Euler-Lagrange equations of our theory to set the 
conditions on 0 and h at the edge. As a result, we are able 
to capture the edge profiles of small membranes whose 
rods are clearly not parallel to the membrane. Presum¬ 
ably, KM theory is amenable to the same approach and 
could thus calculate edge profiles of small membranes. 
KM also view the membrane thickness at the 7 r-wall as a 
boundary condition determined by experiment, whereas 
it is a prediction of our model once physical parameters 
have been set. 

The KM fits to edge and 7 r-wall retardance data re¬ 
ported in Refs. urns (e.g., Fig. 6 of Ref. [42]) are seem¬ 
ingly better than the fits in Figs.Qb) and[ 8 |b) from our 
theory. It should be noted, however, that we use one 
set of parameters to fit data from all membrane radii, 
whereas the KM fits only consider data from a single ra- 
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dius. Our fit to individual profiles are as good as those 
of the KM theory. 


B. Edge ripple fluctuations 

Our free energy Eq. [7] can also be used to investigate 
edge ripple fluctuations of large membranes. First, we 
minimize the free energy over h{x) and 0{pc) to obtain 
the profile for the unperturbed membrane edge. We then 
introduce a small edge ripple u(y) with corresponding 
tangent angle a(y) = d y u(y). We assume that the edge 
profile completely propagates into the membrane inte¬ 
rior, so h(x, y) = h(x — u(y )), and that the rod tilt follows 
the tangent of u(y), so the nematic director changes from 
n(x) = sin h(x)y+cos 6(x)z to n(x, y) = sin a(y) sin 0(x— 
u(y))H + cos a(y) sin Q(x — u(y))y + cos 0{pc — u(y)) z [for 
a schematic of the ansatz, see Fig.[5|a)]. We have to red¬ 
erive the depletion and Frank terms in Eq. [7] to allow for 
gradient terms in the ^/-direction (expression not shown 
here). We expand the ripple tangent angle in Fourier 
components a p : 


a 



(8) 


p is the ripple wavenumber |9]. With the help of a p = 
pu pi where u p s are Fourier components for u(y ), we can 
write the free energy in terms of the small a p s. The free 
energy relative to the state without ripples becomes 

A 77 1 

= 2 E (7[fc, + «[*. e \P 2 ) <4 + 0({«p} 4 ), (9) 

V nr, 


which describes a ID interface with effective line tension 
7 and line bending modulus k. They are given by 


negative at low temperatures, leading to the starfish in¬ 
stability. It is analogous to the chiral line tension term 
in the theory of Langmuir-Blodgett films, which if suffi¬ 
ciently negative, can drive an instability transition from 
a circular film to one with similarly extended arms 
Next we investigate the dynamics of ripple fluctua¬ 
tions. We view the membrane edge as an effective ID 
viscous fluid described by the ripple profile u(y, £), which 
can vary with time. We estimate the Reynolds number of 
this motion to be very small 10 6 -10 4 , so the ripple 
velocity v = d t u obeys overdamped ID hydrodynamics: 

-7710^7 = /ext = /dragH — (13) 

7id is the ID edge viscosity and /dragM is the viscous 
drag force per unit length arising from membrane edge 
motion relative to the bulk solvent [43] . Different mod¬ 
els of membrane-fluid interactions lead to different ex¬ 
pressions for /drag; we see in Appendix [B] that it can 
be largely ignored for ripple wavenumbers probed by our 
experiments. In other words, dissipation of ripple excita¬ 
tions occurs mainly through the membrane rather than 
surrounding solvent since the membrane has much higher 
viscosity. Using = A F/L y — J2 p f p u p for the total 
Hamiltonian density, where A F/L y is given by Eq. [ 9 ] 
and the f p s are an external field formally included to 
calculate the response function, we obtain: 


VidP 2 d t u p = - ( 7 p 2 + up 4 ) u p + f p . 


This leads to the response function 

-1 ( \ 9 f p ( uj ) . 2 . 2 1 4 

= a77) = “ lwr?1DP + + Kp ■ 

The fluctuation-dissipation theorem gives the autocorre¬ 
lation function: 


7 [ft, 0\ = 2 naksT 


f 


dx 


(■ d x h) 2 


\A + (d x h) 2 


K 0 ) 


+ 2 K f dxh [(d x 0) 2 + qd x 0] , ( 10 ) 

Jo 


tv[h 1 6] = 2K f dxh sin 2 0. 
Jo 


( 11 ) 


At thermal equilibrium, the ripple tangent angle compo¬ 
nents take the equipartition values 


Su p u p (^) 


2 k B T 

UJ 


I mXupUp 


2 k B T 1 

ymP 2 u 2 + u p 2 ’ 


where 


7 + Kp 2 

UJ p = 

hlD 


(14) 


is the autocorrelation decay rate. Indeed, temporal ripple 
angle autocorrelations are given by 


( a p) = ~^~2- (!2) 

Note that the term proportional to the chiral twist 
wavenumber q in Eq. [10] is negative for d x 0 < 0. The 
variation of its magnitude with temperature [q(T) is 
temperature-dependent] is the theoretical basis for the 
chiral control of line tension presented in Fig. [5] and 
Ref. |9I. All the other terms are positive-definite, so this 
term must be responsible for the line tension becoming 


(a p (t)a p ( 0 )) = J ^e lu ’ t p 2 S UpUp (uj) = (o£)< 


— 0J v t 


with (a 2 ) in Eq. 
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(15) 


V. EXPERIMENTAL METHODS 

As model rod-like colloids, we use two strains of the fil¬ 
amentous fd bacteriophage: wildtype (wt) and the Y21M 
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mutant [21]. As compared to fd- wt, /d-Y21M has a sin¬ 
gle point mutation in which the 21st amino acid of the 
major coat protein is changed from tyrosine (Y) to me¬ 
thionine (M). Both viruses have the same contour length, 
880 nm, and diameter, 6.6 nm; their persistence lengths 
are 2.8 pm for fd -wt and 9.9 pm for fd- Y21M. They form 
cholesteric phases with opposite handedness: fd -wt forms 
left-handed cholesterics whereas /d-Y21M forms right- 
handed cholesterics. Finally, the chirality of fd -wt is 
temperature-sensitive whereas the chirality of /d-Y21M 
is temperature-independent |9L 

Both viruses are synthesized using standard biological 
protocols [44]. After synthesis, we observe a small por¬ 
tion of viruses that are very long—two and three times 
the nominal length of the virus. We fractionated the 
viruses through the isotropic-nematic phase transition; 
only the isotropic fraction, enriched in nominal-length 
viruses, is kept for this work [9]. These monodisperse 
viruses are then dispersed with concentration c v i ms = 
lmgmL -1 in 20 mM Tris buffer at pH 8.0 and 100 mM 
NaCl. Dextran (500 kDa, Sigma-Aldrich) is used as a 
depletant agent. 

Samples are prepared between glass cover slides and 
coverslips in homemade chambers. A layer of unstretched 
Parafilm is used as a spacer. Slides are coated with poly¬ 
acrylamide brushes to prevent nonspecific binding of the 
viruses with the glass slides and to suppress the deple¬ 
tion interaction between viruses and the glass walls |24| . 
Samples are made airtight using UV-treated glue (Nor¬ 
land Optical). Microscopy observations were performed 
with the inverted microscope Nikon Eclipse Ti equipped 
with an oil immersion objective (1.3 NA, lOOx Plan- 
Fluor). Data is acquired using a cooled CCD camera 
(Andor Clara) for low acquisition rates (below 50 Hz) and 
Phantom v9.1 (Vision Research) for fast acquisition rates 
(above 1000 Hz). 

Sample temperature is tuned between 4 and 60 °C 
with a homemade Peltier module equipped with a 
proportional-integral-derivative temperature controller 
(ILX Lightwave LPT 5910). The temperature-controlling 
side of the Peltier device is attached to a copper ring 
fitted around the microscope objective, which heats or 
cools the sample through the immersion oil. A thermis¬ 
tor, placed in the copper ring adjacent to the sample, en¬ 
abled the proportional-integral-derivative feedback nec¬ 
essary to adjust the temperature. Excess heat is removed 
using a constant flow of room-temperature water. Such 
a device allows us to trigger the starfish instability as 
shown in Fig. [lje). 

The local tilt of the rods with respect to the optical axis 
of the microscope is determined using an LC-Polscope 
(Cambridge Research and Instrumentation) [45]. LC- 
PolScope produces images in which the intensity of each 
pixel is the local retardance D of the membrane. Such 
images can be quantitatively related to the tilting of 
the rods away from the membrane normal (the z-axis 
in Fig. [l]). Rods in the bulk of a membrane are aligned 
along the z- axis, and LC-PolScope images appear black 


in that region. By contrast, for sufficiently large mem¬ 
branes, the bright birefringent ring along the membranes 
periphery indicates local rod tilting as shown in Fig.^a). 
In Fig.|8ja), the LC-PolScope image of a 7r-wall indicates 
that the structure contains twist. 

The time-independent analysis of thermal edge ripple 
fluctuations with DIC optical microscopy yields the line 
tension and the bending rigidity of the edge [9j[46]. The 
acquisition is performed at 1 Hz so that the edge fluctu¬ 
ations are decorrelated. Intensity profile cuts along the 
perpendicular to the edge are fitted by a Gaussian and 
yield the conformation of the edge with subpixel accu¬ 
racy. Each conformation is described in terms of the 
Fourier amplitudes a p (Eq. p|. Averaging over a suffi¬ 
cient number of uncorrelated images gives a fluctuation 
spectrum as shown in Fig. |5|b), where the mean-square 
amplitude (ct p ) is plotted as a function of the wavenum¬ 
ber p. The dynamical analysis of thermal edge ripple 
fluctuations with DIC optical microscopy yields the auto¬ 
correlation decay timescale. The acquisition is performed 
at 3000 Hz. The autocorrelation decay timescale 1 /uj p at 
a given wavenumber p is obtained by fitting the temporal 
autocorrelation function of the Fourier amplitudes by a 
simple exponential (Eq. [I5| ) . Measurements over a suffi¬ 
ciently long time give 1/uj p as a function of p as shown 
in Fig.^c). 

Colloidal membranes can be manipulated using optical 
tweezers. The laser tweezers setup is built around an in¬ 
verted Nikon TE-2000 microscope. A 1064 nm laser beam 
(Compass 1064, Coherent) is projected onto the back fo¬ 
cal plane of an oil-immersion objective (Plan Fluor lOOx, 
NA = 1.3) and subsequently focused onto the imaging 
plane. Using custom Lab VIEW software, multiple trap 
locations were specified and used to stretch and manip¬ 
ulate membranes. Above 2 W of laser power, one can rip 
off smaller membranes from a larger membrane to pro¬ 
duce membranes between 0.5 to 5 pm in diameter. This 
technique is used to study small membranes as shown in 

Fig. [lp). 


VI. DISCUSSION 

The microscopic building components required for as¬ 
sembly of colloidal membranes are monodisperse rod- 
like viruses, non-adsorbing dextran polymer, and poly¬ 
electrolytes to screen electrostatic repulsion. Despite 
their relative simplicity, these building blocks can as¬ 
semble into a myriad of complex structures. Our theory 
demonstrates how their rich properties can emerge from 
hard rods and depletants through three simple entropic 
considerations: depletant excluded volume, rod fluctu¬ 
ations perpendicular to the membrane, and rod twist¬ 
ing as described by the Frank free energy. For example, 
the curved membrane edge with chiral rods arises from 
the competition between depletion, which prefers a cir¬ 
cular vertical edge profile, and Frank elasticity, which 
prefers an exponential edge twisting profile (Fig. [3] and 
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Appendix 0. If depletion is strong enough compared to 
the Frank contribution, achiral virus mixtures will also 
form twisted membranes through spontaneous symmetry 
breaking. Furthermore, our theory predicts that smaller 
membranes, with less distance over which rods can twist, 
have more squared-off edge profiles; this prediction was 
verified by additional experimental data (Fig. [§. De- 
creasing the temperature increases the preferred twist 
wavenumber and consequently decreases the energy of 
the membrane edge, where the twist is greatest. Thus, 
ripple fluctuations, which lengthen the membrane edge, 
increase in amplitude (Fig. [ 5 ]). Eventually, at low enough 
temperatures, edges are energetically preferred and rip¬ 
ples are stabilized in a twisted starfish configuration 
(Fig. |6|. Besides explaining the properties of the mem¬ 
brane edge, our theoretical model can also explain the 
structure of tt- walls. Along membrane edges, a high 
depletion concentration strongly suppresses rod fluctu¬ 
ations perpendicular to the membrane (Fig. |3|. Along 
the middle of 7r-walls, however, large rod fluctuations, 
which can be interpreted experimentally as rod stacking, 
decreases the depletants’ excluded volume and are thus 
favored (Fig. [ 7 ]). This stack of rods with finite thick¬ 
ness physically connects the two partially coalesced mem¬ 
branes and, through depletion, keeps them together. 

All variables in our theory have direct physical mean¬ 
ing. We directly manipulate two of these parameters— 
temperature and depletant concentration—and measure 
several independent physical properties—membrane re- 
tardance (Figs. [4] and [8| and edge fluctuation spectra 
(Fig. §. Theoretical calculations of these properties 
demonstrate respectable agreement with experimental 
measurements while using physically reasonable param¬ 
eter values (Table 0. We use values for the hard-sphere 
depletant size a and fd virus birefringence An that are 
within ^25% of the reported values. We require q(T) to 
have its measured square-root behavior. The Frank elas¬ 
tic constant K and preferred twist wavenumber q(T) are 
^5 times larger than the values measured from viruses 
dispersed in a bulk cholesteric phase without any deple¬ 
tant. However, K and q depend strongly upon the virus 
concentration [23]; membranes condensed by depletants 
have a higher virus concentration than cholesteric sus¬ 
pensions do and thus should have higher K and q. 

Our theory uses a number of assumptions and sim¬ 
plifications. We study the membrane in the continuum 
limit with only two coarse-grained degrees of freedom. 
We ignore rod-rod interactions other than those implicit 
in the phenomenological Frank free energy, whose moduli 
are assumed to be equal and temperature-independent. 
Rod fluctuations perpendicular to the membrane do not 
directly increase the membrane’s volume in the simple 
manner assumed, and while these fluctuations are most 
important at large rod angles 0 ~ 7r/2, their energetic 
cost (Eq. [ 5 ]) was calculated in the small rod angle, small 
fluctuation amplitude limit. In addition, the retardance 
formula was derived for a material of constant thickness 
and optical axis, which does not apply to our membranes. 


We assume a simple ripple ansatz to calculate edge fluc¬ 
tuation spectra, but the actual ripples may have a differ¬ 
ent configuration with lower energy. Yet, despite all of 
these approximations, our model can match experimen¬ 
tal results with quantitative accuracy, indicating that it 
still has value in describing and elucidating properties of 
colloidal membranes. 

The role of depletion and other hard-core interactions 
in colloidal systems has been vigorously investigated from 
many perspectives. Direct excluded volume minimiza¬ 
tion was used to study deplet ion-driven helix formation 
in elastic tubes [47] . Effective entropic potentials be¬ 
tween two anisotropic colloidal particles have been calcu¬ 
lated in depth [481 and have been used to explain various 
self-assembly processes EB EMI]. Free-volume theory 
and theories based on pair distribution functions have 
probed the depletion-induced phase separation of col¬ 
loidal species and have provided relatively sophisticated 
expressions for effective interfacial tensions [ 52 ? ? , [53]. 
However, to our knowledge, the depletion interaction has 
never appeared before as an effective surface tension of 
magnitude naksT explicity. Our system admits this ex¬ 
pression because there is near-complete phase separation 
between the colloids and the depletants and because de¬ 
pletion is strong enough to fix the membrane volume in 
the continuum limit. Our naksT surface tension can be 
related to scaling arguments near the coexistence line in 
Flory-Huggins-de Gennes theory, which proposes an in¬ 
terfacial tension proportional to ksT //f nter , where /inter 
is the thickness of the interface between colloid-rich and 
colloid-poor phases [52j? ? ]. Taking this thickness 

approximately to be the equilibrium rod height fluctu¬ 
ation amplitude bo calculated in our theory, our surface 
tension expression agrees with that obtained by scaling: 
k B T/lf nter ~ k B T/bl naksT. Moreover, the ability 
of our model to quantitatively match and predict exper¬ 
imental results supports the validity of our expression, 
which may guide the design of other colloidal systems 
whose surface tension can be easily tuned by changing 
depletant concentration, depletant size, or temperature. 
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Appendix A: Spontaneous chiral symmetry breaking 
at membrane edges 

As discussed in Section IIV Al and demonstrated in 
Fig. I c), rod height fluctuations are strongly suppressed 
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in membrane edge configurations. We can simplify the 
free energy (Eq. [7]) by enforcing h = tcosO + bo and 
approximating bo = 0. The free energy can then be ex¬ 
pressed in terms of 0 only. In a dimensionless form with 
x = x/t, L x = L x /t, q = g£, and F = F/2natkBTL y , 


F = 




1 + sin 2 0 (d x 6) 2 — cos 0 ) + cos 0(0) 


+ ^ f dx cos 0 (d x 6) 2 — kq sin 6(0). (Al) 

2 Jo 


To investigate the onset of twist, we expand this free 
energy for small 6. To third order, the first integral of 
the Euler-Lagrange equation gives 


y/kdze = -e+ 

zdk 

This equation at x = 0 can be combined with the varia¬ 
tional boundary condition 


kds:0(O ) = -kq - 0(0) + q0 2 ( 0) + T4M 3 (0) 


to obtain 0(0). We first consider q = 0, so F has chiral 
symmetry. We find a twist solution when k < k c = 1, 
where 


0(0) « ±y-(i - &) (A2) 

close to the critical point. When k > 1, only the trivial 
0 = 0 solution exists. If we allow a small nonzero q to 
break the chiral symmetry, a twist solution appears above 

k c . 


m K (as) 

We can integrate the Euler-Lagrange equation to leading 
order and obtain 

6(x) « 0(0) exp(— x/Vk). (A4) 

y/kt acts as a twist penetration depth /twist in analogy to 
smectic phases. Free energy calculations confirm that the 
twist solutions are favored whenever they exist. Thus, 
when q = 0, the phase transition at the k c = 1 criti¬ 
cal point is second-order and spontaneously breaks chiral 
symmetry. Above & c , there is a critical second-order line 
at q = 0. 

We also investigate the edge profile when k <<L It 
is more convenient to write the free energy (Eq. |A1| ) in 
terms of h = h/t = cos 0: 

F = * dx (^iFFF) 2 -h)+ M0) 

+ \ X - kq\[l - M(0). (A5) 


We chose the sign of the square-root in the last term 
assuming 0 > 0, so this expression applies for q > 0. If 
q < 0, then 0 < 0 configurations have lower energy and 
we should choose the opposite sign. The first integral of 
the Euler-Lagrange equation gives 

, 1 * “<^> 2 

A+iw 2 i-i > 2 

This equation at x = 0 can be combined the variational 
boundary condition 

0 = M0) _ 1 

^/i + (^Mo)) 2 

0)d$h(0) _ fe(0) 

1- ^ 2 (°) q JFFFo) 

to obtain a twist solution as a power series in k: 

m^sjFk-^-kq. (A6) 

Solving the Euler-Lagrange equation with k = 0 yields a 
circular profile 


h(x) 



x < 1 — 


x > 1 — 


EM 

2 

fr 2 (o) 

2 


( A7 ) 

However, since cos 0(0) = h( 0) <C 1, the coupling in F ro d 
may be weak. Calculations using the full free energy 
should be performed to check if h = tcosO is a valid 
assumption. 


Appendix B: Estimation of solvent drag during 
ripple fluctuations 


Here we estimate the dissipative forces exerted by the 
solvent on the membrane as it undergoes ripple fluctua¬ 
tions. We approximate the membrane as an infinite 2D 
fluid and apply the analysis of [54], who consider the 
drag force exerted by a subfluid of depth d below the 
fluid plane. For a velocity field v = rx with wavevector 
P = Py, the drag per unit area is #drag = — ? 7 sPcoth (pd) v, 
where r] s is the subfluid viscosity. We estimate a solvent 
depth d ~ 0.3 pm under the membrane where the poly¬ 
mer brush lies. The fluid above the membrane plane 
exerts much less drag because coth pd is a monotonically 
decreasing function of d, so we ignore it. Assuming that 
an effective width l x ~ /twist ~ t of the membrane edge 
moves during the ripple fluctuations, the drag force per 
unit length is approximately / drag = f d xg drag ~ gdrag/*- 
This force modifies the fluctuation autocorrelation decay 
constant (Eq. 14) to 


idr> — 


7 + hip 2 


Vs lx coth (pd)/p + rj ID ’ 
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With rj s « 3mPas from H3 and t]id ~ 300mPaspm 2 
from this work, this change would increase the calcu¬ 
lated values of l/u p in Fig. [5jc) at low wavenumbers 
P < 0.3 pm -1 , but it would not significantly modify our 
fit value for 77 id- For example, at p = 0.3pm -1 , l/u) p 


would be increased ^20%. Moreover, since the measured 
values of l/u) p do not show any increase at small p, this 
analysis may overestimate the drag force, a claim whose 
verification would require a much more sophisticated the¬ 
ory that better captures the ripple geometry and motion. 
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